Boundary Conditions of the Hydro-Cascade Model and Relativistic Kinetic Equations 

for Finite Domains 
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A detailed analysis of the coupled relativistic kinetic equations for two domains separated by 
a hypersurface having both space- and time-like parts is presented. Integrating the derived set of 
transport equations, we obtain the correct system of the hydro-f cascade equations to model the 
relativistic nuclear collision process. Remarkably, the conservation laws on the boundary between 
domains conserve separately both the incoming and outgoing components of energy, momentum and 
baryonic charge. Thus, the relativistic kinetic theory generates twice the number of conservation 
laws compared to traditional hydrodynamics. Our analysis shows that these boundary conditions 
between domains, the three flux discontinuity, can be satisfied only by a special superposition of 
two cut-off distribution functions for the "out" domain. All these results are applied to the case of 
the phase transition between quark gluon plasma and hadronic matter. The possible consequences 
for an improved hydro-l-cascade description of the relativistic nuclear coUisions are discussed. The 
unique properties of the three flux discontinuity and their effect on the space-time evolution of the 
transverse expansion are also analyzed. The possible modifications of both transversal radii from 
pion correlations generated by a correct hydro-fcascade approach are discussed. 
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I. INTRODUCTION 

The modern history of relativistic hydrodynamics 
started more than fifty years ago when L. D. Landau 
suggested [1] its use to describe the expansion of the 
strongly interacting matter that is formed in high energy 
hadronic collisions. Since that time there arose a funda- 
mental problem of relativistic hydrodynamics known as 
the freeze-out problem. In other words, one has to know 
how to stop solving the hydrodynamical equations and 
convert the matter into free streaming particles. There 
were several ways suggested to handle it, but only re- 
cently a new approach to solve the freeze-out problem 
in relativistic hydrodynamics has been invented by Bass 
and Dumitru (BD model) [2] and further developed by 
Teaney, Lauret and Shuryak (TLS model) [3] . These hy- 
dro + cascade models assume that the nucleus-nucleus 
collisions proceed in three stages: hydrodynamic expan- 
sion (hydro) of the quark gluon plasma (QGP), phase 
transition from the QGP to the hadron gas (HQ) and the 
stage of hadronic rescattering and resonance decays (cas- 
cade). The switch from hydro to cascade modeling takes 
place at the boundary between the mixed and hadronic 
phases. The spectrum of hadrons leaving this hypersur- 
face of the QGP-HG transition is taken as input for the 
cascade. 

This approach incorporates the best features of both 
the hydrodynamical and cascade descriptions. It allows 
for, on one hand, the calculation of the phase transition 
between the quark gluon plasma and hadron gas using 
hydrodynamics and, on the other hand, the freeze-out 



of hadron spectra using the cascade description. This 
approach allows one to overcome the usual difhculty 
of transport models in modeling phase transition phe- 
nomenon. For this reason, this approach has been rather 
successful in explaining a variety of collective phenomena 
that has been observed at the CERN Super Proton Col- 
lider (SPS) and Brookhaven Relativistic Heavy Ion Col- 
lider (RHIC) energies. However, both the BD and TLS 
models face some fundamental difRcultics which cannot 
be ignored (see a detailed discussion in [4]). Thus, within 
the BD approach the initial distribution for the cascade 
is found using the Cooper- Frye formula [5] , which takes 
into account particles with all possible velocities, whereas 
in the TLS model the initial cascade distribution is given 
by the CMi-Oj(f formula [6,7], which accounts for only those 
particles that can leave the phase boundary. As shown 
in Ref. [4] the Cooper-Frye formula leads to causal and 
mathematical problems in the present version of the BD 
model because the QGP-HG phase boundary inevitably 
has time-like parts. On the other hand, the TLS model 
does not conserve energy, momentum and number of 
charges and this, as will be demonstrated later, is due 
to the fact that the equations of motion used in [3] are 
incomplete and, hence, should be modified. 

These difficulties are likely in part responsible for the 
fact that the existing hydro-|-cascade models, like the 
more simplified ones, fail to explain the HBT puzzle [8], 
i.e. the fact that the experimental HBT radii at RHIC 
are very similar to those found at SPS, even though the 
centre of mass energy is larger by an order of magni- 
tude. Therefore, it turns out that the hydro-fcascade 
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approach successfully parameterizes the one-particle mo- 
mentum spectra and their moments, but does not de- 
scribe the space-time picture of the nuclear collision as 
probed by two-particle interfcrometry. 

The main difficulty of the hydro -I- cascade approach 
looks similar to the traditional problem of freczo-out in 
relativistic hydrodynamics [6,7]. In both cases the do- 
mains (subsystems) have time-like boundaries through 
which the exchange of particles occurs and this fact 
should be taken into account. In relativistic hydrody- 
namics this problem was solved by the constraints which 
appear on the freeze-out hypersurface and provide the 
global energy- momentum and charge conservation [6,7,9]. 
A generalization of the usual Boltzmann equation which 
accounts for the exchange of particles on the time-like 
boundary between domains in the relativistic kinetic the- 
ory was given recently in Rcf. [4] . It was shown that the 
kinetic equations describing the exchange of particles on 
the time-like boundary between subsystems should neces- 
sarily contain the (5-likc source terms. From these kinetic 
equations the correct system of hydro-|-cascade equations 
to model the relativistic nuclear collision process was de- 
rived without specifying the properties of the separating 
hypersurface. However, both an explicit switch off cri- 
terion from the hydro equation to the cascade one and 
the boundary conditions between them were not consid- 
ered in [4] . The present work is devoted to the analysis of 
the boundary conditions for the system of hydro-|-cascade 
equations. This is necessary to formulate the numerical 
algorithm for solving the hydro-|-cascade equations. 

The paper is organized as follows. In Sect. 2 a brief 
derivation of the set of kinetic equations is given and 
source terms are obtained. In Sect. 3 the analog of the 
collision integrals is discussed and a fully covariant for- 
mulation of the system of coupled kinetic equations is 
found. The relation between the system obtained and 
the relativistic Boltzmann equation is also considered. 
The correct equations of motion for the hydro -|- cas- 
cade approach and their boundary conditions are ana- 
lyzed in Sect. 4. There it is also shown that the existence 
of strong discontinuities across the space-like boundary, 
the time-like shocks, is in contradiction with the basic 
assumptions of a transport approach. The solutions of 
boundary conditions between the hydro and cascade do- 
mains for a single degree of freedom and for many degrees 
of freedom are discussed in Sect. 5 and 6, respectively. 
The conclusions are given in Sect. 7. 



the hypersurface S*. The distribution fimction (f)in{x,p) 
for t < t*{x) is assumed to belong to the "in" domain, 
whereas (j)outix,p) denotes the distribution function of 
the "out" domain for t > t*{x){see Fig.l). In this work 
it is assumed that the initial conditions for 0i„(x,p) are 
given, whereas on S* the function 0ouiix,p) is allowed 
to differ from (j)in{x,p) and this will modify the kinetic 
equations for both functions. For simplicity we consider 
a classical gas of point-like Boltzmann particles. 
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FIG. 1. Schematic two dimensional picture of the bound- 
ary hypersurface E* (solid curve). Arrows show the exter- 
nal normal vectors. The light cone NOP is shown by the 
dash-dotted line. The point F divides E* into the time-like 
(OF) and space-like [FP) parts. 

Similar to Ref. [10] we derive the kinetic equations for 
(j)in{x,p) and (po'uiiX'P) from the requirement of parti- 
cle number conservation. Therefore, the particles leaving 
one domain and crossing the hypersurface S* should be 
subtracted from the corresponding distribution function 
and added to the other. Now consider the closed hyper- 
surface of the "in" domain, Ax^ (shown as the contour 
ABDE in Fig.l), which consists of two semi-planes an 
and (7(2 of constant time tl and t2, respectively, that are 
connected from t\ to t2 > tl by the arc BD of the bound- 
ary Tj*{tl,t2) in Fig.l. The original number of particles 
on the hypersurface at\ is given by the standard expres- 
sion [10] 



II. DRIFT TERM FOR SEMI-INFINITE DOMAIN 

Let us consider two semi-infinite domains, "in" and 

"out", separated by the hypersurface S* which, for the 
purpose of presenting the idea, we assume to be given in 
(3-1-1) dimensions by a single valued function t = t*{x) = 
Xq{x). The latter is assumed to be a unique solu- 
tion of the equation J^*{t,x) = (a switch off criterion) 
which has a positive time derivative doT*{t*,x) > on 



Nr = - j d^^^ p^ ct>in{x,p), (1) 



where rfE^ is the external normal vector to an and, 
hence, the product p'^dT,^ < is non-positive. It is clear 
that these particles can cross either hypersurface 0-42 or 
T,*{tl,t2). The corresponding numbers of particles are 



2 



as follows 

N2 



NT = 

loss 



A* n ^ 



i{x,p), (2) 



e(p''dE,). (3) 



E»(tl,t2) 

The 0-function in the loss term (3) is very important 
because it accounts for the particles leaving the "in" do- 
main (see also discussion in [6,9]). For the space- like 
parts of the hypersurface S*(tl,t2) which are defined 
by negative sign ds"^ < of the squared line element, 
ds^ = dt* {x)^—dx'^ , the product p'^dT,^ > is always pos- 
itive and, therefore, particles with all possible momenta 
can leave the "in" domain through the T,*{tl,t2). For 
the time-like parts of S*(il,t2) (with sign ds'^ > 0) the 
product p'^dY^u can have either sign, and the 0-function 
cuts off those particles which return to the "in" domain. 

Similar one has to consider the particles coming to the 
"in" domain from outside. This is possible through the 
time-like parts of the hypersurface T,*{tl,t2), if the par- 
ticle momentum satisfies the inequality —p"dT,i, > 0. In 
terms of the external normal (iS^ with respect to the 
"in" domain (this normal vector is shown as an arrow 
on the arc BD in Fig.l and will be used hereafter for all 
integrals over the hypersurface T,*{tl,t2)) the number of 
gained particles 

d^p 



N* ■ = - 

gain 



dT, 



r <Pout{x,p) ei-p^di:,) (4) 



E*(tl,f2) 

is, evidently, non-negative. Since the total number of 
particles is conserved, i.e. N2 = Ni - A^;*sg + 
one can use the Gauss theorem to rewrite the obtained 
integral over the closed hypersurface A.t'' as an integral 
over the 4-volume Ax'* (area inside the contour ABDE 
in Fig.l) surrounded by Aa;^ 

J d^x^ pf^ 4>in{x,p) = fdj:^^p^ 



d'^x^ pf" (l)in{x,p) = ^ 
Ax" E*(tl,t2) 

[^in{x,p) - (l)out{x,p)]Q{-p''di:^) . (5) 



P" 



Note that in contrast to the usual case [10], i.e. in the 
absence of a boundary S*, the right-hand side (rhs) of 
Eq. (5) does not vanish identically. 

The rhs of Eq. (5) can be transformed further to a 4- 
volume integral in the following sequence of steps. First 
we express the integration element dJ^^ via the normal 
vector n* as follows {dx^ > 0, for j = 1, 2, 3) 

d'^n = n*dx^dx'^dx^; n* = 6^0 - ^^^(1 - S^lo) , (6) 
where S^i, denotes the Kronecker symbol. Then, using 

t2 



where the 4-dimcnsional volume is a direct prod- 
uct of the 3- and 1-dimensional volumes I]*(tl,t2) and 
{t2 — tl), respectively. Evidently, the Dirac (5- function 
allows us to extend integration in (7) to the unified 4- 
volume V^ — Ax^ U of Ax* and (the volume 
is shown as the area ABCE in Fig.l). Finally, with the 
help of notations 



eo„t = e(t-r(x)); ei„ = i-eo„t 



(8) 



it is possible to extend the left hand side (Ihs) integral in 
Eq. (5) from Ax^ to V^. Collecting all the above results, 
from Eq. (5) one obtains 



d^x^ Gin p^ <i>in = J d^x^ p^n; X 



d^p 



[^in - <^o„t] Qi-P^K) 6{t - t* (X)) . (9) 

Since the volumes Ax'* and V^ are arbitrary, one obtains 
the kinetic equation for the distribution function of the 
"in" domain 

©m dn (j)in{x,p) = Cin{x,p) + 

p^n;[<l>in{x,p) - cl>out{x,p)] ei-p^nt) S{t - t*{x)) . (10) 

Note that the general solution of Eq. (9) contains an 

arbitrary fimction Cj„(x,p) (the first term in the rhs of 
(10)) which identically vanishes while being integrated 
over the invariant momentum measure d^p/po. Such a 
property is typical for a collision integral [10], and we 
shall discuss its derivation in the subsequent section. 

Similar one can obtain the equation for the distribution 
function of the "out" domain 

Oout P^ cj)out{x,p) = Cout{x,p) + 

p^n;[0,„(x,p) - 0™t(x,p)] e(pX) 5{t - t*{x)) , (11) 

where the normal vector n* is given by (6). Note the 
asymmetry between the rhs of Eqs. (10) and (11): for 
the space-like parts of hypersurface S* the source term 
with 0(— p^n*) vanishes identically because p'^n* > 0. 
This refiects the causal properties of the equations above: 
propagation of particles faster than light is forbidden, 
and hence no particle can (re)enter the "in" domain. 



III. COLLISION TERM FOR SEMI-INFINITE 
DOMAIN. 

Since in the general case 0i„(x,p) ^ (f>ouL{x,p) on E*, 
the (5-like terms in the rhs of Eqs. (10) and (11) cannot 
vanish simultaneously on this hypersurface. Therefore, 



the identity Jdt5{t - t3) = 1 for the Dirac (5-function functions 6*, = ein\i:- and e*„t = Goutl^' ^ 



tl 



with tl <t3 < t2, we rewrite the rhs integral in (5) as 



J 

S«(tl,t2) 



d^x 6{t-t*{x)) n* 



(7) 



do not vanish simultaneously on S* as well. Since there 
is no preference between "in" and "out" domains it is 
assumed that 



e* = 6* 



(12) 



3 



but the final results are independent of this choice. 

Now the collision terms for Eqs. (10) and (11) can 
be readily obtained. Adopting the usual assumptions 
for the distribution functions [11.10,12], one can repeat 
the standard derivation of the collision terms [10,12] and 
get the desired expressions. We shall not recapitulate 
this standard part, but only discuss how to modify the 
derivation for our purpose. First of all, one has to start 
the derivation in the Ax^ volume of the "in" domain and 
then extend it to the unified 4- volume = Ax'^ U 
similarly to the preceding section. Then the first part of 
the collision term for Eq. (10) reads 

CUx,p) = el ilG[ct>in,<t>in] - lL[<t>in, <f>in]) , (13) 

Jg[0a, ^b] = IJ D^P <Pa{p) 4>b{Pi) W^ppJp'p; , (14) 

Il[<^a, M = \j D^P Mp) Mpi) W^PP.bX , (15) 

where the invariant measure of integration is denoted by 
D^P = and W„„ uy,y is the transition rate 

in the elementary reaction with energy-momentum con- 
servation given in the form p'^ + pf^ = p''^ + p^ . The 
rhs of (13) contains the square of the Oi„-function be- 
cause the additional 0i„ accounts for the fact that on 
the boundary hypersurface S* one has to take only one 
half of the traditional collision term (due to Eq. (12) only 
one half of S* belongs to the "in" domain). It is easy to 
understand that on S* the second part of the collision 
term (according to Eq. (12)) is defined by the collisions 
between particles of "in" and "out" domains 

Cll{x,p) = 6i„eo„t {lG{<i>inAout\ - lL[(t>in, '/'out]) • (16) 

Combining (10), (13) and (16), one gets the kinetic equa- 
tion for the "in" domain 

Q^n d^(j),nix,p) = C/„(a;,p) + Cl^{x,p)+p>'nl x 
[(/.„(a;,p) - (j>out{x,p)] e(-p'^0 5{t - t*{x)) . (17) 

The kinetic equation for the "out" domain can be de- 
rived similarly and then it can be represented in the form 

©out P^ d^(f>out{x,p) = Cl^t{x,p) + C'„^t{x,p) +p^n*^ X 

[(f>in{x,p) - (t>out{x,p)] eip^K) S{t - t*{x)) , (18) 

where the evident notations for the collision terms 
CLt = {lG[(f>out,(f>out] - lL[(f'out,(f>out]) and = 

6m0o«t (-^g[0 out') 4^in out 1 4^in 

]) are used. 

The equations (17) and (18) can be represented also in 
a covariant form with the help of the function J^*{t,x). 
Indeed, applying the definition of the derivative of the 
implicit function to d^t*{x), one can rewrite the exter- 
nal normal vector (6) as n* = dij,J^* {t,x)/ doJ^* {t,x). 
Now using the inequality doT*{t*,x) > and the fol- 
lowing identities S{T*{t,x)) = S{t - t* (x)) / doT* {t* , x) , 
Qa = Q{Sa J^*{t,x)) one can write Eqs. (17) and (18) 
in a fully covariant way 



cf^Aix.p) = &a{x,p) + c'Jix,p) + p^d^r X 

\'\>in{x,p) - <^out{x,p)\ ^{Sa p'^d^T*) 5{T*{t,x)) , (19) 

where the notations A € m. Sin — {A £ out, 
= 1) are introduced for "in" ("out") domain. 
For the continuous distribution functions on E*, i.e. 
0oMt|s» = 4'in\Y.' , the source terms on rhs of Eqs. (17) 
and (18) vanish and one recovers the Boltzmann equa- 
tions. Moreover, with the help of the evident relations 

e,„ = Qout = S{J'*{t, x)) 5^^*(i, x) , (20) 
+ + Cl^, + Cli, = 7g[$, - II $] , (21) 

where ^{x,p) = 6j„ 4>in{x,p) + Qout (t>out{x,p), one can 
get the following result summing up Eq. (17) and (18) 

9^ $(a:,p) =/g[$,$] -/l [$,$]. (22) 

In other words, the usual Boltzmann equation follows 
from the system (19) automatically without any assump- 
tion about the behavior of (pin and (pout on the boundary 
hypersurface S*. Also Eq. (22) is valid not only under 
condition (12), but for any choice < < 1 obeying 
Eq. (8). 

In fact the system (19) generalizes the relativistic ki- 
netic equation to the case of the strong temporal and 
spatial inhomogeneity, i.e., for (j)in{x,p) ^ 4>out{x,p) on 
S*. Of course, one has to be extremely careful while dis- 
cussing the strong temporal inhomogeneity (or disconti- 
nuity on the space-like parts of S*) such as the so called 
time-like shocks [13,14] because, as shown in the subse- 
quent section, their existence contradicts the usual as- 
sumptions [11,10,12] adopted for distribution functions. 

From the system (19) it is possible to derive the macro- 
scopic equations of motion for the energy-momentum 
tensor by multiplying the corresponding equation with 
p" and integrating it over the invariant measure. Thus, 
Eq. (19) generates the following expression {T'^ = 

/$^pWa(x,p)) 

^" = 1^ p''C'J{x,p) + P^'p^d^T* X 

{<Pin{x,p) - <l>out{x,p)] Q{Sa p'dpT*) S{r'{t,x)). (23) 

Similar to the usual Boltzmann equation the momentum 
integral of the collision term C/„ vanishes due to its sym- 
metries [10], but it can be shown that the integral of 
the second collision term C// does not vanish because it 
involves two different distribution functions. 

The corresponding system of equations for the con- 
served current N'^ = J ^ p^^(pA{x,p) can be obtained 
by direct integration of the system (19) with the invariant 
measure 

eAd^N^ = J ^ J^d^r'[,i>in{t,x) - ct>out{t,x)] X 

Q{SAP'dpJ'*)H^*it,x)). (24) 
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The above equation does not contain the contribution 

from antiparticles (just for simplicity), but the latter can 
be easily recovered. Note that in contrast to (23) the 
momentum integral of both collision terms vanish in Eq. 
(24) due to symmetries. 



IV. CONSERVATION LAWS AT E* 

It is clear that Eqs. (19), (23) and (24) remain valid 
both for finite domains and for a multiple valued function 
t = t*{x) as well. To derive the whole system of these 
equations in the latter case, one has to divide the function 
t*{x) into the single valued parts, but this discussion is 
beyond the scope of this paper. Using Eqs. (19), (23) and 
(24) we are ready to analyze the boundary conditions on 
the hypersurfacc S*. The simplest way to get the bound- 
ary conditions is to integrate Eqs. (23) and (24). Indeed, 
integrating (23) over the 4- volume V^j, (shown as the area 
ABCD in Fig. 2) containing part pT,* of the hypersur- 
face S* , one obtains the energy- momentum conservation. 
Before applying the Gauss theorem to the Ihs of (23), we 
note that the corresponding ©^-function reduces the 4- 
volumc VpY: to its part which belongs to the ^-domain. 
The latter is shown as area ALMD [BCML) for A&in 
{A e out) in Fig. 2. Then in the limit of a vanishing 
maximal distance A ^ between the hypersurfaces AD 
and BC in Fig. 2, the volume integral of the Ihs of Eq. 
(23) can be rewritten as the two integrals / ddfiT'^ : 
the first integral is performed over the hypersurfacc 
shown as an arc LM in Fig. 2, and the second integral 
reduces to the same hypersurfacc but taken in the op- 
posite direction, i.e. the ML arc in Fig. 2. Thus, the 
volume integral of the Ihs of Eq. (23) vanishes in this 
limit, and we obtain 

= Jd^x QaO^ {n''{x,p)) ^Jd^^^ p''Ci'{x,p) + 
J rf4a;^5(^*(i,x))pVa^J-* X 



[(pinix.p) - 4>outix,p)] 0(5.4 p''dpT*) 



(25) 



Similarly to Sect. II, in the limit A ^ the second inte- 
gral on the rhs of (25) can be reexpressed as an integral 
over the closed hypersurfacc. Since the latter is arbi- 
trary, then Eq. (25) can be satisfied, if and only if the 
energy-momentum conservation occurs for every point of 
the hypersurfacc S* 

Tt:^ d,J^*{t*,x)=TZ± d,J^*{t*,x), (26) 
T'^l^ j'^ P^p''cI>a{x,p) e(± pPdpJ^*) . 

In deriving (26) from (25) we used the fact that the 4- 
volume integral of the second collision term vanishes 



for finite values of distribution functions because of the 
Kronecker symbols. The results for the conserved cur- 
rent follows similarly from Eq. (24) after integrating it 
over the 4- volume V^j^ and taking the limit A ^ 

Kn± 9,r' {t* , x) = N^^ d^r' it* , x) , (27) 

N^^^ J'^ p^<Pa{x.p) e(± ppdpj^*) . 

The fundamental difference between the conservation 
laws (26), (27) and the ones of usual hydrodynamics is 
that the systems (26) and (27) conserve the quantities of 
the outgoing from (Sa = 1) and incoming to (Sa = — 1) 
"in" domain particles separately, whereas in usual hy- 
drodynamics only the sum of these contributions is con- 
served. 
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FIG. 2. Schematic two dimensional picture of the integra- 
tion contour to derive the boundary conditions (25) - (27) 
between the "in" and "out" domains. In the limit of a van- 
ishing maximal distance A — » between the hypersurfaces 

AD and BC, both of these hypersurfaces arc reduced to the 
part pE* (an arc LM) of the boundary E* between domains. 

The trivial solution of Eqs. (26) and (27) corresponds 
to a continuous transition between "in" and "out" do- 
mains 



4>out{x,p) (l)in{x,p) 



(28) 



This choice corresponds to the BD model [2] . The BD 
model gives a correct result for an oversimplified kinetics 
considered here. However, in the case of the first or- 
der phase transition (or a strong cross-over) which was a 
prime target of the hydro-|-cascade models [2,3] the situ- 
ation is different. In the latter case the speed of soimd ei- 
ther vanishes (or becomes very small) [15,16] and, hence, 
the rarefaction shock waves become possible [17-19]. The 
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reason why the rarefaction shocks may exist Hcs in the 
anomalous thermodynamic properties [19] of the media 
near the phase transition region. In other words, on the 
boundary between the mixed and hadronic phases the 
rarefaction shocks are mechanicahy stable [19], whereas 
the compression shocks are mechanically unstable. This 
is also valid for the vicinity of the generalized mixed 
phase of a strong cross-over. 

One important consequence of the shock mechanical 
stability criterion is that the stable shocks necessarily are 
supersonic in the media where they propagate. The latter 
means that the continuous rarefaction flow in the region 
of phase transition is mechanically unstable as well, since 
a rarefaction shock, if it appears, propagates inside the 
fluid faster than the sound wave and, hence, it should 
change the fluid's state. Due to this reason the unstable 
hydrodynamic solutions simply do not appear [20,21]. 

Applying these arguments to the BD model, one con- 
cludes: for the first order phase transition or strong cross- 
over the sound wave in the (generalized) mixed phase 
may be unstable and the strong discontinuities of the 
thermodynamic quantities are possible [17-19]. The lat- 
ter corresponds to the non-trivial solution of the conser- 
vation laws (26) and (27), which allows a discontinuity of 
the distribution function on two sides of the hypersurface 
S* . Since there is twice the number of conservation laws 
compared to the usual hydrodynamics, it is impossible, 
as shown below, to build up the nontrivial solution of 
Eqs. (26) and (27), if the distribution functions on both 
sides of the hypersurface E*, i.e. (jiin and cf)out^ are taken 
to be the equilibrimn ones. 

Consider first the space-like parts of the hypersurface 
E*. Then Eqs. (26) and (27) for Sa = -I vanish 
identically because of the inequality p^5pJF*(f*, .t) > 0, 
whereas for Sa = ^ Eqs. (26) and (27) recover the usual 
hydro dynamical conservation laws at the discontinuity. 
However, it can be shown that the existence of strong dis- 
continuities across the space-like hypersurfaces, the time- 
like shocks [13,14], is rather problematic because it leads 
to a contradiction of the basic assumptions adopted for 
the distribution function, even though the conservation 
laws (26) and (27) are formally fulfilled. 

Indeed, according to the Bogolyubov's classification 
[1 1] , a one-particle treatment can be established for a typ- 
ical time At which, on one hand, should be much larger 
than the collision time tcoIU and, on the other hand, it 
should be much smaller than the relaxation time TReiax 

TColl «At« (29) 

Similar to the usual Boltzmann equation (see also dis- 
cussions in [11,12]), in deriving the collision terms of 
Eq. (19) we implicitly adopted the requirement that the 
distribution function does not change substantially for 
times At less than the relaxation time TReiax- However, 
at the discontinuities on the space- like parts of E*, sug- 
gested in [13,14], the distribution function changes sud- 
denly, i.e. At = 0, and the left inequality (29) cannot 



be fulfilled at the time-like shock. Therefore, according 
to the Bogolyubov's classification [11], such a process, 
which is shorter than the typical collision time, belongs 
to a prekinetic or chaotic stage and, hence, cannot be 
studied at the level of a one-particle distribution func- 
tion. It would instead require the analysis of a hierar- 
chy of A''-particle distribution functions, where A'' is the 
number of particles in the system. Thus, the existence 
of time-like shocks contradicts the adopted assumptions 
for a one-particle distribution. Their existence should be 
demonstrated first within the higher order distributions. 
This statement applies to several papers published by the 
Bergen group during the last few years where time-like 
shocks were attenuated in time using a phenomenological 
quasi-kinetic approach [22] . For the same reason, the use 
of equilibrium values for temperature and chemical po- 
tential in an attenuated time-shock is rather problematic 
for time scales shorter than rcoii- Note, however, that the 
discontinuities at the time-like parts of E* (usual shocks) 
have no such restrictions and, hence, in what follows we 
shall analyze only these discontinuities. 

V. BOUNDARY CONDITIONS AT E' FOR A 
SINGLE DEGREE OF FREEDOM 

Now we have to find out whether it is principally pos- 
sible to obtain the nontrivial solution of systems (26) 
and (27) using the parts of equilibrium distributions on 
the time-like segments of the hypersurface E*. To sim- 
plify the presentation, first we consider the same kind of 
particles in both domains. As usual in relativistic hy- 
drodynamics it is convenient to transform the coordinate 
system {t*{x); x) into the special local frame (t^(xi); x^) 
being, for definitencss, the rest-frame of discontinuity be- 
tween the distributions (/)i„ and (j)out and defined as fol- 
lows: the X-axis should coincide with the local external 
normal vector to the hypersurface E*, whereas y- and 
z-axes should belong to the tangent hypcrplane to E*. 
In this case the external normal vector to the time-like 
parts of E* has a simplest form: n* = ( 0; 0; 0), 

and one can readily check that the value of the deriva- 
tive diJ^^ plays an important role in the conservation 
laws (26) and (27) only through the the cut-off © func- 
tion. Then, like in the theory of usual relativistic shocks 
[20,19,21], it can be shown that equations for y- and z- 
components of system (26) degenerate to the identities 
because of the symmetries of the energy-momentum ten- 
sor. Therefore, the number of independent equations at 
the discontinuity is 7: a switch off criterion and six inde- 
pendent equations out of systems (26) and (27) {t- and 
x-equations (26) and one equation (27) for two choices of 
Sa = {-1;+1}). 

On the other hand the number of unknowns is 6 only: 
temperature T^^ and baryonic chemical potential of 
the "in" domain, temperature T*^^ and baryonic chemical 
potential /Lt*„j of the "out" domain, the collective veloc- 
ity w*„ of the "in" domain particles, and the collective 
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velocity of the particles of "out" domain, which, evi- 
dently, should be coUinear to the normal vector n* in the 
rest-frame of the discontinuity. Thus, a formal counting 
of equations and unknown shows that it is impossible to 
satisfy the conservation laws (26) and (27), if the distri- 
bution functions on both sides are the equilibrium ones. 

The last result means that instead of a traditional dis- 
continuity we have to search for a principally new bound- 
ary conditions on the hypersurface S*. The analysis 
shows that there are two of such possibilities with the 
equilibrium distribution function in the "in" domain and 
a special superposition of two cwi-o/f equilibrium distri- 
butions for the "out" domain. The first possibility is to 
choose ^o«t as follows: 



(30) 



i.e. the distribution of outgoing particles from the "in" 
domain (the first term in the rhs of Eq. (30)) is con- 
tinuous on the hypersurface S*, whereas the distribu- 
tion of the incoming to "in" domain particles (the sec- 
ond term in the rhs of Eq. (30)) has a discontinuity 
on S* which, however, conserves the energy, momentum 
and baryonic charge because of the following boundary 
conditions {v = {0; 1}) 

^in— i'-^in: f^iji- ^out— {'-^outT ^out^ ^oui^ • (31) 
^in- C^in^ f^in^'^in) ^ ^out- {"^ouf l^out^^out) ■ (32) 



The above choice of boundary conditions at S* allows 
to reduce the systems (26) and (27) for = 1 to the 
identities, and, hence, from the systems (26) and (27) 
there remain only three independent equations (31), (32) 
for 5*^ = — 1. Alone with a switch off criterion, these four 
equations can be solved now for six independent vari- 
ables, and, consequently, the two variables can be chosen 
free for a moment. Thus, we showed that both the outgo- 
ing and incoming parts of the distribution function (30) 
can be chosen as the equilibrium ones, but with different 
temperatures, chemical potentials and non-zero relative 
velocity vl^i = {vl^^ - <J/(1 - «*„*<„) with respect to 
the distribution function 

Note a principal difference between this discontinuity 
and all ones known before: the "out" domain state con- 
sists, in general, of two different subsystems (fluxes) that 
have their individual hydrodynamic parameters. It is 
clear that it is impossible to reduce three of those hy- 
drodynamical parameters of one flux to those three of 
other flux because there arc only two free variables out 
of six. Thus, together with the "in" domain flux there 
are in total three fluxes involved in this discontinuity, 
and; therefore, it is appropriate to name it a three flux 
discontinuity in order to distinguish it from the ordinary 
shocks that are defined maximum by two fluxes. 

The outgoing component of the distribution (30) co- 
incides with the choice of the boundary conditions sug- 
gested in the TLS model [3], whereas the equations (31) 



and (32) are missing in this model. This is the reason 
why the TLS model suffers from the energy, momentum 
and charge non-conservation. It is also necessary to note 
that the lower values of the temperature T*^^ < T*^ and 
baryonic chemical potential < y^*„, which are typ- 
ical for the rarefaction process considered in [3], should 
be compensated by an extra flow from the incoming par- 
ticles to the "in" domain, i.e. v*^i should be opposite to 
the external normal vector n* in the rest-frame of the 
three flux discontinuity. Therefore, such a discontinuity 
is analogous to the compression shock wave in relativis- 
tic hydrodynamics, and, hence, it cannot appear in the 
rarefaction process for any of the hadronic species con- 
sidered in Ref. [3]. 

Similarly, one can find another non-trivial solution of 
the systems (26) and (27) which corresponds to opposite 
choice to Eq. (30) 



0m ( T*„, mL, <J ei-p'd,Ti) + 



(33) 



i.e., the incoming to the "in" domain component of the 
distribution above (the first term in the rhs of Eq. (33)) 
is continuous on hypersurface S* , but the outgoing from 
the "in" domain component has a discontinuity on S* 
which obeys the following conservation laws (i^ = {0; 1}): 



rpLl> /rp* * 

in+ in^ Pin 7 ' 
(T* II* - 



'in) 



— T^OUt+ C^OUtT I^OUtT'^OUt) 7 (34) 

— ^out+ (^o«t) Mout' "^otit) • (35) 



It is clear that both the outgoing and incoming compo- 
nents of the distribution (33) can be chosen as the equilib- 
rium distribution functions. A simple analysis of the sys- 
tem (34), (35) shows that for T*^^ < T*^ and fj,*^^ < 
the relative velocity Vrei in the local frame should be 
coUinear to the external normal vector n*, and, hence, 
such a discontinuity is analogous to the rarefaction shock 
wave in the relativistic hydrodynamics. Thus, in contrast 
to the TLS choice, Eq. (33) should be used as the ini- 
tial conditions for the "out" domain while studying the 
rarefaction process of matter with anomalous thermody- 
namic properties. 

Now we are ready to discuss the question how the non- 
trivial solutions (30) and (33) will modify the system 
of the hydro-l-cascade equations (19), (23) and (24). In 
what follows we shall assign the hydrodynamic equations 
to the "in" domain and the cascade ones to the "out" 
domain (the opposite case can be considered similarly). 
Applying Eqs. (30), (31) and (32) to the "in" Eqs. (23) 
and (24) and to the "out" Eq. (19), one obtains the 
following system: 



Qout P" <i>out{x,p) = Ci^t{x,p) + Ciit{x,p) , (38) 



• ^IJ, -'-in 
Hn <Jii ^^in 



(36) 
(37) 
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i.e., due to the boundary conditions (30) - (32) the 6- 
hkc terms have disappeared from the original system of 
equations. Also it is clear that the source term in the rhs 
of Eq. (36) does not play any role because it is finite on 
the hypersurface S* and it vanishes everywhere outside 
of S*. 

In order to obtain the system of hydro+cascade equa- 
tions (36) - (38) for the non-trivial solution defined by 
Eqs. (33) (35), the hydrodynamic description has to 
be extended to the outer e- vicinity (e 0) of the hyper- 
surface S* 



(39) 
(40) 



which in practice means that for Eqs. (33) (35) one has 
to start solving the cascade equation (38) just bit inside 
of the "out" domain in order to get rid of the ^-like term 
in (38) and move this term to the discontinuity on the 
hypersurface S*. 

The remarkable feature of the system of hy- 
dro+cascade equations (36) - (40) is that each equation 
automatically vanishes outside of the domain where it is 
specified. Also, by the construction, it is free of the prin- 
cipal difiiculties of the BD and TLS models discussed 
above. The question how to conjugate the three flux dis- 
continuity with the solution of the hydro equations (36) , 
(37), (39) and (40) will be discussed in the next section. 



VI. BOUNDARY CONDITIONS AT E* FOR 
MANY DEGREES OF FREEDOM 

In order to apply the above results to the description 
of the QGP-HG phase transition that occurs in relativis- 
tic nuclear collisions it is necessary to take into account 
the fact that the real situation differs from the previous 
consideration in two respects. The first one is that in the 
realistic case inside the "in " domain there should ex- 
ist the QGP, whereas it should not appear in the "out" 
domain. Of course, the discussion of the QGP kinetic 
theory is much more complicated problem and it lies 
far beyond the subject of this work. For our purpose 
it is sufficient to generalize the equations of motion (36) 
- (40) inside domains and the conservation laws (26) and 
(27) between these domains to the realistic case. Such 
a generalization can be made because in the case of the 
QGP-HG phase transition there will be also an exchange 
of particles between the "in" and "out" domains which 
has to be accounted by the (5-like source terms in the 
transport equations. The only important difference from 
the formalism developed in the preceding sections is that 
QGP must hadronize while entering the "out" domain, 
whereas the hadrons should melt while entering the "in" 
domain. Note, however, that in relativistic hydrodynam- 
ics one has to assume that all reaction, i.e. the QGP 
hadronization and melting of hadrons in this case, occur 



instantaneously. Under this assumption one can justify 
the validity of the equations of motion (36) - (40) and the 
conservation laws between QGP and HG on the bound- 
ary S*. 

The second important fact to be taken into account 
is that some hadrons have the large scattering cross- 
sections with other particles and some hadrons have the 
small cross-sections, and because of that the hadrons of 
both kinds participate in the collective flow differently. 
A recent effort [23,24] to classify the inverse slopes of 
the hadrons at SPS lab energy 158 GeV-A led to the 
conclusion that the most abundant hadrons like pions, 
kaons, (anti)nucleons, A hyperons e.t.c. participate in 
the hadron rescattering and resonance decay till the very 
late time of expansion, whereas 17 hyperons, J/ip and 
■0' mesons practically do not interact with the hadronic 
media and, hence, the freeze-out of their transverse mo- 
mentum spectra {kinetic freeze-out) may occur just at 
hadronization temperature Th- Therefore, the inverse 
slopes of the fi, J/iJj and ip' particles is a combination of 
the thermal motion and the transversal expansion of the 
media from which those particles are formed. 

These results for the O baryons and (p mesons were ob- 
tained within the BD and TLS models, whereas for the 
J/tjj and V' mesons it was suggested in Refs. [23,24] for 
the first time. Later on these results were refined further 
in Ref. [25] by the simultaneous fit with the only one free 
parameter (the maximal value of transversal velocity) of 
the measured fl [26,27], J/ip and tp' [28] transverse mo- 
mentum spectra in Pb-|-Pb collisions at 158 GeV-A that 
are frozen-out at hadronization temperature Th ■ The ex- 
perimental situation with the (j) mesons at SPS is, unfor- 
tunately, not clarified yet because the results of the NA49 
[29] and NA50 [30] Collaborations are not in agreement. 
The analysis of the transverse momentum spectra of f2 
hyperons [31,32] and (f) mesons [31] reported by the STAR 
Collaboration for energies y/s = 130 A-GcV in Refs. [33] 
and [34], respectively, and for = 200 A-GeV in Ref. 
[32] shows that this picture remains valid for RHIC en- 
ergies as well. 

It is easy to find out that for such particles like (p, 
il, J/ip and ip' which are weakly interacting with other 
hadrons the distribution function (pout should coincide 
with <p,n 



I'out 



(41) 



where, in contrast to (30), there is no incoming compo- 
nent of the distribution because the non-interacting par- 
ticles cannot rescatter and, hence, change their velocity. 
Note also that a small modification of the incoming part 
of J/ ^ momentum distribution due to decay of heavier 
charmonia in the "out" domain can be safely neglected. 
Remarkably, the cascade initial condition (41) exactly co- 
incides with the one used in the TLS model. Therefore, 
the main TLS conclusions [3] on the <p mesons and f2 
hyperons remain unchanged, whereas for hadrons with 
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large scattering cross-sections the TLS conclusions may 
change significantly. 

Omitting the contributions of weakly interacting 
hadrons from the components of the energy-momentum 
tensor and baryonic 4-current, one can generalize the 
boundary conditions (26) and (27) on the hypcrsurface 
S* between the domains, and formulate the energy- 
momentum and charge conservation laws in terms of the 
cut-off distribution functions. For dcfinitcncss wc shall 
consider the first order phase transition between QGP 
and hadronic matter through out the rest of this work. 
The case of second order phase transition can be ana- 
lyzed similarly. Thus, in terms of the local coordinates 
{t}^{xL)',XL), introduced in Sect. 5, the conservation laws 
(26) and (27) can be generalized as follows {f = {0; 1}) 

Q=q,q,... 

(1 - a,) Th±{ TL, Zh ■ mL, O = 

H=7V,K,... 

E Th''±{T^^t,ZH-t^tt,vfut), (42) 

H=v,K,... 

Q=q,q,... 

(1-"?) E ^H±{Tin, Zh ■ H*„, v*J = 

H=TT,K,... 

E Nh±{Ttut,ZH-f^tut,vtut) (43) 

H=v,K,... 

where a<j is the volume fraction of the QGP in a mixed 
phase, the Q-sums of the energy-momentum tensor and 
baryonic 4-current components, denoted as 

T^±^l^ P'P'' M^^P) ©(± P'dp^n - (44) 
N^±^l^ P" ZqcI>q{x,p) e(± p'^d.T*) , (45) 

run over all corresponding degrees of freedom of QGP. 
Similarly, the _ff-sums run over all hadronic degrees of 
freedom. In Eqs. (42) and (43) Zq and Zh denote the 
baryonic charge of the corresponding particle species. 

Now from Eqs. (44) and (45) it is clearly seen that the 
correct hydro-)-cascade approach requires the knowledge 
of a more detailed information on the microscopic proper- 
ties of QGP than it is usually provided by the traditional 
equations of state. To proceed further we, however, shall 
assume that those components are known. The general 
approach to calculate the angular and momentum inte- 
grals in Eqs. (44) and (45) was developed in Ref. [7] and 
was applied to the massive Boltzmann gas description in 
[35,7]. 

The important difference between conservation laws 

(42), (43) and (26), (27) is that in the "out" domain 
the temperature T~^^, chemical potential and rela- 
tive velocity of the incoming to S* hadrons should 



differ from the corresponding quantities T^^^, and 
v'^^^ of the outgoing from E* particles, and both sets 
should differ from the quantities T*^, ii*^ and v*^ of the 
"in" domain. In order to prove this statement it is nec- 
essary to compare the number of equations and number 
of unknowns for the two distinct cases, namely, (i) if the 
initial state is in the mixed QGP - HG phase, and (ii) if 
the initial state belongs to the QGP. 

In case (i) there are 10 equations and 10 unknowns: 

• The equations are as follows: 6 conservation laws 
from Eqs. (42) and (43); value of the initial energy; value 
of the initial baryonic density; the relation between initial 
temperature T^^ and the baryonic chemical potential 
taken at the phase boundary; and the switch off criterion. 

• The unknowns arc as follows: three temperatures 
^out> ^ou*; three chemical potentials /i^^j, ji^^^.; 

three velocities Vg^^, defined in the rest frame of 
a discontinuity; and the QGP fraction volume Oq. Thus, 
in this case one can find a desired solution of the system 
of ten transcendental equations, which is the most gen- 
eral form of the three flux discontinuity introduced by 
Eqs. (30 ) - (32). 

To complete the solution of hydro equations (36), (37), 
(39) and (40) one has to find out the value of velocity 
w*„ from the system of ten transcendental equations dis- 
cussed above. Then this velocity defines an ordinary dif- 
ferential equation dx\/dt*j^ = —v^^ for the hypersurface 
S* in the rest frame of matter of the "in" domain, which 
has to be solved simultaneously with the hydro equations. 

If initial state belongs to the interior of the QGP phase, 
case (ii), then the usual hydro solution will be valid till 
the systcan reaches the boundary with the mixed phase, 
from which the non-trivial discontinuity described by 
Eqs. (42) and (43) will start on. Now it is clear what 
are the distinctive features from the previously consid- 
ered case: in contrast to case (i) on the boundary with 
the mixed phase the volume fraction of QGP is fixed to 
unit Uq = 1\ the energy and baryonic charge densities are 
not independent anymore, but are completely defined by 
the temperature and baryonic chemical potential, which, 
in addition, are connected by the entropy conservation 
for the continuous hydro solution in QGP. 

Therefore, in case (ii) there are 9 equations and 9 un- 
knowns, which are as follows: 

• The equations are: 6 conservation laws from Eqs. 
(42) and (43); temperature dependence of the baryonic 
chemical potential /Lt|„ = l^i^n^in) ^'^^ *° entropy con- 
servation; the relation connecting temperature T*^ and 
baryonic chemical potential , since they belong to the 
phase boundary; and the switch off criterion. 

• On the other hand the unknowns, except for the 
fixed volume fraction Uq = 1, are the same as in case (i). 
Thus, again the number of unknowns matches the num- 
ber of equations, and the procedure to solve the system 
of hydro equations (36), (37), (39) and (40) simultane- 
ously with the boundary conditions (42) and (43) is the 
same as in case (i). 
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Now it is appropriate to discuss the switch off criterion 
J-*{t, a;) = in more details. By the construction of the 
hydro+cascade approach, the cascade treatment should 
be applied since the very moment, where the hydrody- 
namics starts to lose its applicability: according to the 
original assumption the hydro equations (36), (37), (39) 
and (40) work well inside of the 4- volume surrounded by 
the hypersurface E* and in the outer e-vicinity (e 0) of 
E* [see also a discussion after Eq. (38)], whereas just out- 
side of this domain the thermal equilibrium dismantles 
and, hence, one has no right to use the cwi-oj^' equilibrium 
distributions interior of the "out" domain. Consequently, 
a switch off criterion should be formulated solely for some 
quantity defined in the outer e-vicinity of hypersurface 
S*, and it has to define the bounds of applicability of 
thermal equilibration and/or hydrodynamic description. 
Note that in the BD and TLS models this did not matter 
because both groups kept the cascade initial conditions 
as close as possible to the output of hydro. However, in 
the case of the three flux discontinuity on the time-like 
parts of hypersurface E* the proper use of the switch off 
criterion plays a decisive role in the construction of the 
mathematically correct hydro-|-cascade solution (see also 
a discussion of the freeze-out criterion in Refs. [6,7]). It is 
clear that, in contrast to the BD and TLS formulations, 
the switch off criterion may generate a very sizable effect 
while applied to interior of hadronic phase. It is so, since 
even a small difference (just a few MeV) between the 
temperature T*^, which belongs to the phase transition 
region, and temperatures T~^^ and Tj^j of the "out" do- 
main may lead to a tremendous flow of outgoing hadrons 
because of the enormous latent heat of the QGP. 



VII. CONCLUDING REMARKS 

In the preceding sections we have derived the system 
of relativistic kinetic equations which describe the parti- 
cle exchange between two domains separated by the hy- 
persurface of arbitrary properties. We showed that the 
usual Boltzmann equation for the following sum of two 
distributions '^{x,p) = Q^n 4>in{x,p) + 9o«t (l)oui{x,p) 
automatically follows from the derived system, but not 
vice versa. Integrating the kinetic equations we derived 
the system of the hydro+cascade equations for a single 
degree of freedom. Remarkably, the conservation laws 
on the boundary between two domains conserve the in- 
coming and outgoing components of the energy, momen- 
tum and baryonic charge separately, and, hence, there 
is twice the number of conservation laws on the separat- 
ing hypersurface compared to the usual relativistic hy- 
drodynamics. Then we showed that for a single degree of 
freedom these boundary conditions between domains can 
be satisfied only by a special superposition of two cut-off 
equilibrium distributions for the "out" domain. Since 



the obtained discontinuity, in contrast to usual shocks 
defined by two fluxes, has three irreducible fluxes, hence, 
it is named a three flux discontinuity. It was also shown 
that the TLS-like choice of the boundary conditions, in 
contrast to expectation of [3], corresponds to an analog 
of the compression shock in traditional hydrodynamics, 
and, therefore, it cannot be used to model the rarefaction 
process. 

Then we showed that existence of the time-like shocks 
[13,14], formally rederived by this formalism, contradicts, 
nevertheless, to the usual assumptions adopted for the 
one-particle distributions and, hence, the solution of this 
problem requires the analysis of higher order distribu- 
tion functions. Therefore, in the rest of the paper we 
concentrated on a detailed analysis of the discontinuities 
at the time-like hypersurfaces, i.e. the space-like shocks 
in terms of Refs. [13,14]. These results were then general- 
ized to a more realistic case, namely, if the mixed QGP- 
HG phase is assigned to the "in" domain and hadrons 
exist in the "out" domain. Such a generalization also re- 
quired the exclusion of the hadrons with the small scat- 
tering cross-section (like fl, J/^p and tp' particles) from 
the boundary conditions between domains. As we showed 
in the preceding section, the presence of the first order 
phase transition makes the resulting system of transcen- 
dental equations more complicated than in the case of a 
single degree of freedom. 

It turns out, that a minimal number of variables in 
this discontinuity is either 9 or 10, depending on the lo- 
cation of the initial state on the phase diagram, and, 
therefore, on the hadronic side the three flux discontinu- 
ity should have two different flows with their own tem- 
peratures, chemical potentials and collective velocities. 
The found solution has a number of unique features in 
comparison with usual shocks: 

• this discontinuity may generate a very strong, ex- 
plosive like, flow of outgoing particles from the "in" 
domain, first, because a huge latent heat of QGP is 
involved, and, second, due to an extra momentum 
associated with the cut-off distribution. Indeed, 
considering the outgoing component of the distri- 
bution (pout Q{p^diT1) for massless pions in the 
frame where this function maximally resembles the 
non-cut Boltzmann distribution, i.e. in the rest- 
frame of the latter, one finds a nonvanishing collec- 
tive velocity = -^i^!^. Here v„ = ^ (|u<t| < 1 
for time- like parts of E*) denotes the transversal 
radius velocity in this frame. 

• the strong explosive flow of outgoing particles is lo- 
calized at the time-like parts of the hypersurface 
E*, whereas at the space-like parts of E* there will 
be a continuous flow. It is even possible that for 
some choice of parameters the space-like boundary 
may be absent. 

• the particle density of outgoing pions will strongly 
depend on the speed of the transversal radius ex- 
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pansion. Thus, for massless pions the particle den- 
sity found according to the Eckart definition [10] is 

Ptt = ^"^4°"'^ Vi^ ~ + i.e. it is smaller 

for all ?v > ^1 than the thermal particle density 
PniT^t)- Therefore, the two particle correlations 
off the low particle density regions should be re- 
duced. Since, the situation v^^ » —1 is typical 
for the beginning of the transversal expansion [3], 
then the main contribution to the transversal pion 
correlations will come from the later times of ex- 
pansion. Thus, it is possible that the space-time 
region which defines the side and out pion correla- 
tion radii will be essentially more localized both in 
space and time than in traditional hydrodynamic 
solutions. 

• since there arc two fluxes in the "out" domain, they 
will interact with each other. The resulting distri- 
bution should be, of course, found by the cascade 
simulations, but it is clear that the fastest of them 
will decelerate and the cold one will reheat. Besides 
the possibility to accelerate or decelerate the out- 
going transversal flow more rapidly than in the BD 
and TLS models, the three flux discontinuity may 
naturally generate some turbulence patterns in the 
"out" domain. 

Taking into account all these features alone with the fact 
that neither the BD nor TLS boundary conditions have 
such a strong discontinuity, we conclude that the three 
flux discontinuity opens a principally new possibility not 
only to resolve the HBT puzzle [8] , but to study some new 
phenomena, like a turbulence pattern, associated with a 
new kind of shock, a three flux discontinuity, in relativis- 
tic hydro-|-cascade approach. 

Despite the reasonably good description of the one- 
particle spectra of the most abundant hadrons, even 
such sophisticated model as the TLS one badly overes- 
timates both of the transverse radii measured by pion 
interferometry like other hydrodynamic models. This is 
a strong indication that the hydro part of all existing hy- 
dro+cascade and hydrodynamic models requires an es- 
sential revision. How this revision will affect the present 
BD and TLS results is unclear at the moment, but the 
solution of the HBT puzzle [8] should serve as a good test 
for the correct picture of the space-time evolution dur- 
ing the post-hadronization stage. The additional tests 
for the correct hydro+cascade equations should be the 
reproduction of three recently established signals of the 
deconfinement phase transition, i.e. the pion kink [36,37] 
seen at lab energy of ^ 30 GeV-A , the K~^/tt~^ peak 
at the same lab energy [36] (the strangeness horn) and 
the plateau [38] in the inverse slope of the transverse 
momentum spectra at the whole range of the SPS ener- 
gies (the step in caloric curves) measured by the NA49 
Collaboration [39,40]. Also it is necessary to check other 
predictions of the Statistical Model of the Early Stage 
[36], namely the anomalies in the entropy to energy fluc- 



tuations [41] (the " fink") and in strangeness to energy 
fluctuations [42] (the "tooth"), because both the "fink" 
and "tooth" may be sensitive to the turbulence behavior 
due to energy dissipation. 

Note, however, that the completion of this task re- 
quires an additional research of the hydro+cascade ap- 
proach. First, it is necessary to develop further the mi- 
croscopic models of the QGP equation of state in order 
to find out the required by Eqs. (42) - (45) components 
of the cut-ofl^ energy-momentum tensor and baryonic 4- 
current. This can be done, for example, within the phe- 
nomenological extensions [43^5] of the Hagedorn model. 
Second, a similar problem for hadrons should be solved 
as well, otherwise, as we discussed in preceding section, 
the switch off criterion from the hydro to cascade can- 
not be formulated correctly within the hydro-|-cascade 
approach. And, finally, for practical modeling it is nec- 
essary to formulate a mathematical algorithm to solve 
simultaneously the system of hydro-|-cascade equations 
(36) - (40) with the boundary conditions (42) and (43) 
between the hydro and cascade domains. These prob- 
lems, however, should be considered elsewhere. 

Acknowledgments. The author is thankful to 

D. Blaschkc, W. Gassing and L. W. Phair for valuable 
comments, and to P. Huovinen and B. R. Schlei for stim- 
ulating discussions. Also I appreciate very interesting 
discussions on this subject with L. P. Csernai, E. V. 
Shuryak and D. Teaney. The author thanks the Institute 
for Nuclear Theory at the University of Washington for 
its warm hospitality and the Department of Energy for 
the partial support during the completion of this work. 
Also the partial financial support of the DFG grant No. 
436 UKR 17/13/03 is greatly acknowledged. 



[1] L. D. Landau, Izv. Akad. Nauk Ser. Fiz., 51 (1953). 

[2] S. Bass and A. Dumitm, Phys. Rev. C61, 064909 (2000); 

[3] D. Teaney, J. Lauret and E. V. Shuryak, Phys. Rev. Lett. 
86, 4783 (2001); and nucl-th/0110037 (2001). 

[4] K. A. Bugacv, Phys. Rev. Lett. 90, 252301 (2003). 

[5] F. Cooper and G. Frye, Phys. Rev. DIO, 186 (1974). 

[6] K. A. Bugaev, Nucl. Phys. A606, 559 (1996). 

[7] K. A. Bugaev and M. I. Gorenstein, nucl-th/9903072. 

[8] M. Gyulassy, Lect. Notes Phys. 583, 37 (2002). 

[9] K. A. Bugaev, M. I. Gorenstein and W. Greiner, J. of 
Phys. G25, 2147 (1999); Heavy Ion Phys. 10, 333 (1999). 
[10] S. R. de Groot, W. A. van Leeuwen and Ch. G. van 
Weert, Relativistic Kinetic Theory (North-Holland Pub- 
lishing, Amsterdam 1980). 
[11] N. N. Bogolyubov, Problems of Dynamic Theory in Sta- 
tistical Physics (Gostekhizdat, Moscow 1946) (in Rus- 
sian). [Reprinted in Studies in Statistical Mechanics, 



11 



vol. 1 (J. de Boer and G. E. Uhlenbeck, eds., North- 
Holland Publishing, Amsterdam 1962).] 

[12] R. Balescu, Eqmltbrium and Nonequilibrium Statistical 
Mechanics (Kriger Publishing, Malabar, Florida 1991). 

[13] L. P. Csernai, Zh. Eksp. Teor. Fiz. (Russ.) 92, 379 (1987). 

[14] L. P. Csernai, Sov. Phys. JETP 65, 216 (1987). 

[15] C. M. Hung, E. V. Shuryak, Phys. Rev. Lett. 75, 4003 
(1995). 

[16] C. M. Hung, E. V. Shuryak, Phys. Rev. C57, 1891 (1998) 
and references therein. 

[17] L. van Hove, Z. Phys. C21, 93 (1983). 

[18] M. Gyulassy et. al, Nucl. Phys. B237, 477 (1984). 

[19] K. A. Bugaev, M. I. Gorenstein and V. I. Zhdanov, Z. 
Phys. C39, 365 (1988); K. A. Bugaev and M. I. Goren- 
stein, J. of Phys. G13, 1231 (1986) and Z. Phys. C43, 
261 (1989). 

[20] L. D. Landau and E. M. Lifshitz, Fluid Mechanics (Perga- 

mon. New York, 1979); Hydrodynamics (Nauka, Moscow, 

1986) (in Russian). 
[21] for more details and references see K. A. Bugaev et al., 

Phys. Rev. D40 (1989) 2903. 
[22] the whole list of references can be found in L. P. Csernai 

et al, hep-ph/0401005. 
[23] K. A. Bugaev, M. Gazdzicki and M. L Gorenstein, Phys. 

Lett. B523, 255 (2001). 
[24] K. A. Bugaev, .1. of Pliys. G28, 1981 (2002). 
[25] M. I. Gorenstein, K. A. Bugaev and M. Gazdzicki, Phys. 

Rev. Lett. 88, 132301 (2002). 
[26] E. Anderson et al. (WA97 CoUab.), Phys. Lett. B433, 

209 (1998). 

[27] F. Antinori et al. (WA97 CoUab.), J. of Phys. G27, 375 
(2001). 

[28] M. C. Abreu et al. (NA50 Collab.), Phys. Lett. B499, 85 
(2001). 



[29] S. V. Afanasev et al. (NA49 CoUab.), Phys. Lett. B491, 

59 (2000). 

[30] M. C. Abreu et al. (NA50 CoUab.), Nucl. Phys. A661, 
534c (1999). 

[31] K. A. Bugaev, M. Gazdzicki, M.L Gorenstein, Phys. Lett. 

B544, 127 (2002); and Phys. Rev. C69, 017901 (2003). 
[32] C. Suire et al. (STAR CoUab.), Nucl. Phys. A715, 470 

(2003). 

[33] G. van Buren (STAR CoUab.), talk given at QM2002 
(2002); and J. Adams et al. (STAR CoUab.), nucl- 
ex/0307024. 

[34] C. Adler et al, (STAR CoUab.), Phys. Rev. C65, 

041901(R) (2002). 
[35] L. P. Csernai, Z. Lazar and D. Molnar, Heavy Ion Phys. 

5, 467 (1997). 

[36] M. Gazdzicki and M. L Gorenstein, Acta Phys. Polon. 

30, 2705 (1999). 
[37] M. Gazdzicki, Z. Phys. C66, 659 (1995) and J. of Phys. 

G23, 1881 (1997). 
[38] M.I. Gorenstein, M. Gazdzicki and K.A. Bugaev, Phys. 

Lett. B567, 175 (2003). 
[39] S. V. Afanasiev et al (NA49 CoUab.), Phys. Rev. C66 

054902 (2002). 

[40] C. Ah et al (NA49 CoUab.), J. Phys. G30, S119 (2004). 
[41] M. Gazdzicki, M. I. Gorenstein and St. Mrowczynski, 

hep-ph/0304052. 
[42] M.I. Gorenstein, M. Gazdzicki and O. S. Zozulya, hep- 

ph/0309142. 

[43] M.I. Gorenstein, V. K. Petrov and G. M. Zinovjev, Phys. 

Lett. B106, 327 (1981). 
[44] J. Lotossier, J. Rafelski and A. Tounsi, Phys. Lett. B328, 

499 (1994). 

[45] D. B. Blaschke and K. A. Bugaev, nucI-th/0311021. 



12 



